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Superconducting electrical circuits can be used to study the physics of cavity quantum electrody- 
namics (QED) in new regimes, therefore realizing circuit QED. For quantum information processing 
and quantum optics, an interesting regime of circuit QED is the dispersive regime, where the de- 
tuning between the qubit transition frequency and the resonator frequency is much larger than 
the interaction strength. In this paper, we investigate how non-linear corrections to the disper- 
sive regime affect the measurement process. We find that in the presence of pure qubit dephasing, 
photon population of the resonator used for the measurement of the qubit act as an effective heat 
bath, inducing incoherent relaxation and excitation of the qubit. Measurement thus induces both 
dephasing and mixing of the qubit, something that can reduce the quantum non-demolition aspect 
of the readout. Using quantum trajectory theory, we show that this heat bath can induce quantum 
jumps in the qubit state and reduce the achievable signal-to-noise ratio of a homodyne measurement 
of the voltage. 

PACS numbers: 03.65.Yz, 42.50.Pq, 42.50.Lc, 74.50.+r, 03.65.Ta 



I. INTRODUCTION 

Cavity quantum electrodynamics (QED) is a unique 
tool to study the interaction between light and matter 
at its most fundamental level [l], Q ■ Most interesting is 
the regime of strong coupling where the frequency asso- 
ciated with light-matter interaction is greater than all 
relaxation rates [!, Q. For example, in this regime, ex- 
periments with very high-Q cavities have been able to 
resolve quantum jumps and time-resolved collapse of the 
cavity field @,Q. 

With the strong coupling regime easily accessible, su- 
perconducting electrical circuits offer distinctive advan- 
tages to the study of light-matter interaction 0, H, [0, 
I 111. [l2l [HI, [IH , something which has been realized exper- 
imentally with charge [ll| , flux [l|| E3] , and phase |l8l | 
superconducting qubits. Although the present work ap- 
plies to all physical realizations of cavity or circuit QED, 
here we will focus on superconducting charge qubits cou- 



pled to a transmission line resonator [lj, [15[ . Because 
the qubit can be very strongly coupled to the transmis- 
sion line in this system, it opens the possibility to study 
new regimes of cavity QED. For example, the strong dis- 
persive regime was theoretically studied in Ref. [l9[ and 
experimentally investigated in Ref. [2(| ■ 

In circuit QED, readout of the qubit is done by irra- 
diating the resonator with photons at, or close to, the 
bare resonator frequency while the qubit is strongly de- 
tuned from the resonator. Information about the state of 
the qubit is then encoded in the phase and amplitude of 
the field transmitted and reflected from the resonator. 
In principle, increasing the amplitude of the measure- 
ment drive, hence the photon population of the resonator, 
should increase the rate at which information is gained 



about the qubit. For example, in Ref. [14|, l2l| it was 
estimated that by filling the resonator with the critical 
photon number n cr i t = A 2 /4p 2 , where A is the frequency 
detuning between the qubit and the resonator, and g 
is their interaction strength, one would reach signal-to- 
noise (SNR) ratios of ~ 200, even taking into account 
realistic amplifier noise. Such SNR would easily lead to 
single-shot readout in this system [22| . However high 
SNR have not been experimentally observed. 

In previous these work, the conclusions for the SNR 
were obtained by analyzing the qubit-resonator Hamilto- 
nian in the dispersive approximation, which is a perturba- 
tive expansion of the Jayncs-Cummings Hamiltonian to 
second order in g/ A. However, this approximation fails 
as the number of photon in the resonator increases. As 
n cr it is approached one should expect higher order terms 
in the perturbative expansion, with corresponding non- 
linearities, to be important . While the dispersive approxi- 
mation has been shown to be very accurate in understand- 
ing experimental results for circuit QED at moderate pho- 
ton number population [U, 01, HI M HI HE HI Ha , it 
should break down as the photon population is increased. 
Understanding these corrections is important if we are 
to gain more insights in the measurement process and its 
effect on the qubit. 

Here, this is done by pus hing the dispersive approxi- 
mation used in Ref. [lj, \2M to higher order. These re- 
sults not only apply to circuit QED but also to cavity 
QED and more generally to any physical situation where 
a two-level system is dressed by an oscillator. Examples 
arc quatronium [H, HO] or flux qubits [10, HH, HH cou- 
pled to bifurcating oscillators for readout purposes. In 
particular, the authors of Ref. [HJ find that the qubit 
relaxation rate is strongly enhanced when the non-linear 
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oscillator is in its high-amplitude state compared to its 
low-amplitude state, results which are at least qualita- 
tively consistent with those presented here. 

In Sec. IIIII we find a unitary transformation that ex- 
actly diagonalizes the Jaynes-Cummings Hamiltonian. 
Expanding this to higher orders in g/A allows us to de- 
rive results which are valid for for higher photon numbers 
(but still less then n C rit.)- This transformation is then 
applied on the Hamiltonians describing coupling of the 
resonator and qubit to environmental degrees of freedom. 
Taking advantage of the large separation in energy scale 
in the dispersive regime for the relevant qubit and res- 
onator bath frequency, we obtain in Sec. lIVI a Markovian 
Lindblad-type master equation for the system that takes 
into account higher order dispersive corrections and po- 
tential frequency- variations in the environmental spectral 
densities. 

One of the most important results of this paper is ob- 
tained in Sec. [V] There, building on Ref. |2l| . we elimi- 
nate the resonator degree of freedom from the resonator- 
qubit master equation found in Sec. IIVI and derive an 
effective master equation for the qubit 's reduced density 
operator. This is possible when n < {n K , n cr ;t.} where n K 
is a new maximum photon number set by the cavity de- 
cay rate over the strength of the non-linearity. This effec- 
tive master equation contains the measurement-induced 
dephasing found in Refs. [l9|, [2l|, [23| and the novel ef- 
fect of photon-number dependent qubit relaxation and 
dephasing rates. For example we show how finite photon 
population of the resonator acts as an effective heat bath 
on the qubit. This effective model is shown numerically 
to be very accurate in reproducing the dynamics of the 
full Jaynes-Cummings model. By comparing results ob- 
tained with the linear dispersive approximation, we find 
that non-linear effects become important even for photon 
occupation number significantly below n, cr j t . In Sec. IVI1 
the photon-dependent qubit mixing and dephasing rates 
are discussed in more details and the analytical results 
are compared to numerical calculations. 

In Sec. I VIII a quantum trajectory equation describing 
the evolution of the qubit and resonator under homodyne 
measurement is obtained, as well as a reduced qubit quan- 
tum trajectory equation. This is used to investigate the 
measurement, where we show that the achievable SNR is 
decreased substantially by the non-linear effects. 



II. CIRCUIT QED 

In circuit QED, a superconducting charge qubit is fab- 
ricated inside a transmission line resonator. This system 
is illustrated in Fig. [TJ Focussing on a single mode of the 
resonator, the system Hamiltonian describing this circuit 
takes the Jaynes-Cummings form [l4| 

H S = H + hgl+ (2.1) 




FIG. 1: (Color online) Schematic layout and lumped element 
version of the circuit QED implementation. A superconduct- 
ing charge qubit (green) is fabricated inside a superconducting 
ID transmission line resonator (blue). 

where we have defined 

Hp, = hu r a'a + hu a — (2-2) 

I± = a'a~ ± acr+. (2-3) 

In this expression, uj r is the frequency of the mode of in- 
terest of the resonator, u) a the qubit transition frequency 
and g the qubit-resonator coupling. The operators a^' 
and a± are the creation and annihilation operators for 
the photon field and the qubit. 

Logical operations and readout of the qubit can be 
achieved by applying a microwave signal on the input 
port of the resonator. Choosing a frequency that is close 
to the resonator frequency tu r corresponds to a readout 
of the qubit's state, while frequencies that are close to u a 
can be used to control the qubit [3, EH, HH • This can be 
modeled by the Hamiltonian 

H d = Y J ^{^{t)a^e~ l ^ t + el^ae 1 ^ 1 ) . (2.4) 

k 

Taking into account that signals of different amplitude, 
frequency and phase can be sent simultaneously to the 
input port of the resonator, efc(i) is the amplitude of the 
k th drive and ujk its frequency. In this paper, we will be 
more particularly interested in looking at the effect of a 
measurement (fc = to) drive and will have only this drive. 

A. Damping 

The effect of coupling to environmental degrees of free- 
dom can be described by the master equation [33| 

g=^[H,g}+KV[a}g + ll V[<j-}g+^V[a z }g, (2.5) 

where H is the total Hamiltonian of the system including 
drives 

H = H S + H d (2.6) 

and V[L]g = (2Lgtf - L^Lg - gL^L) /2. In the above 
expression, k is the resonator rate of photon loss, 71 the 
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qubit energy decay rate and 7 V is the qubit rate of pure 
dcphasing. 

This master equation is obtained in the Markov approx- 
imation which assumes that the spectral density of the 
environment is frequency- independent. For high-quality 
factor systems, like high-Q transmission line resonators 
or (most) superconducting qubits, this approximation is 
accurate as the system is probing the environment in a 
very small frequency bandwidth. As we will see in Sec. lIVl 
when going to the dispersive approximation, it can be im- 
portant to take into account the frequency dependence of 
the environment. 



III. DISPERSIVE EFFECTS ON THE 
HAMILTONIAN 

In the limit that dctunning between the cavity and 
the qubit is large, no energy is exchange. In this situa- 
tion, the interaction is said to be dispersive. In analyz- 
ing this interaction, it is convenient to diagonalize the 
Jaynes-Cummings Hamiltonian Eq. (|2.ip by using a uni- 
tary transformation. In this transformed frame, the new 
effective qubit and photon operators are combinations of 
the bare qubit and photon operators. In this sense, the 
qubit acquires "a photon part" and vice versa. This leads 
to, for example, the Purcell effect where a qubit can de- 
cay through the photon decay channel [13, HH ■ 



A. Dispersive Jaynes-Cummings Hamiltonian: 
Linear regime 

In the limit where |A| = \uj a — uj r \ 3> g, the Jaynes- 
Cummings Hamiltonian (|2.1[) can be approximately di- 
agonalized using the unitary transformation 



D 



Linear 



(3.1) 



with A = g/A a small parameter. Using the relation 



e- xx He xx = H + \[H,X]+'— [[H, X] , X] + ■ ■ • (3.2) 

to second order in A, it is simple to obtain the effective 
Hamiltonian describing the dispersive regime 

Linear 



A 2 



Tno r a* 



2g\ 



aa+ 2 



y+G(A 2 ). 

(3.3) 



The qubit transition frequency is shifted by a quantity 
proportional to the photon population 2g\( x a^a). Alter- 
natively, this shift can be seen as a qubit dependent pull 
of the resonator frequency ui r — » u r ± gX. 

As a result, shinning microwaves at the input port 
of the resonator at a frequency close to uj r and mea- 
suring the transmitted signal using standard homodync 
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FIG. 2: (Color online) Trace distance between the simula- 
tion of the full master equation and the dispersive linear 
model presented in Ref [211 ] versus time, a) The initial state 
is |0}(|g) + \e))/V2 and g/2ir = 50MHz, A/2n = 2GHz, 
n crit . = 400, k/2tt = 2.5MHz, 7i/2tt = 0.1MHz, 7^/2tt = 
0.3MHz. A measurement drive of amplitude e m /2ir = 10MHz 
of shape 0.5e m (tanh[(i — fj,l)/a] + 1) with 2n/j,i = 10.0/is 
and 2na = 10.0/xs is applied, corresponding to a mean of 
(a^a) ~ 34 ~ 0.08n cr it. photons, b) Maximum of the trace 
distance vs measurement amplitude. The vertical line indi- 
cates the measurement amplitude used in panel a). 



techniques serves as measurement of the qubit [14j, Il3 . 
l2ll I23I |2~4| . In this approximation, this corresponds to 
a quantum non-demolition (QND) measurement of the 
qubit [ll, [El. 

For this measurement scheme, increasing the number 
of photons in the input beam increases the intensity of the 
output signal which would overcome the noise introduced 
by the amplifier. However the Hamiltonian (|3.3p is only 
valid for a mean photon populations n <C n CT it. = 1/4A 2 
(as is shown by a Taylor expansion of the exact result 
Eq. (|3.8|) ] which means we need to consider non-linear 
corrections to Eq. (|3.3p to understand the dynamics as n 
approaches n CT it.- 

To see the breakdown of the linear approximation, we 
have numerically calculated the time-dependent evolu- 
tion of the system master equation under measurement 
(see Fig. [2] for details) both in the linear dispersive ap- 
proximation, using the approach described in Ref. [21] , 
and with the full non-dispersive Jaynes-Cummings model. 
To compare these results we plot the trace distance 



(for{Px,Q2) = ^Tr (\px - Tr r (g 2 ) |) 



/ £ ((^-(^ 

i£{x,y,z} 



(3.4) 



where p\ is the reduced qubit density matrix found us- 
ing the linear dispersive model presented in Ref [2lT ] , and 
Tr r (£2) is the trace over the resonator of the total density 
matrix of the system found by simulation of the complete 
master equation (|2.5[) . This trace distance is the geomet- 
rical distance between two Bloch vectors, and ranges from 
to 2, with when the two qubit states are the same and 
2 when they are opposite on the Bloch sphere. 

We plot, in figure [2] a), the trace distance for a mea- 
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surcmcnt amplitude which is slowly turned on to reach an 
amplitude corresponding to n w 0.08?i cr it. photons in the 
resonator. Even for this small number of measurement 
photons (compared to ?Vrit.) the trace distance is non- 
negligible which implies the breakdown of the dispersive 
approximation. Moreover, we plot in figure Eb) the max- 
imum of the trace distance (over the simulation time), 
as a function of the maximum measurement amplitude 
e m . Clearly, the trace distance gets worse as the ampli- 
tude is increased. The importance of this effect depends 
on the various parameters entering in the simulation but 
the results shown here are typical. It is clear from these 
numerical results, it is important to take into account 
higher order terms in the dispersive approximation. 



B. Dispersive Jaynes-Cummings Hamiltonian: 
Exact transformation 



Following the derivation presented in Appendix 1X1 
which is similar in spirit to the approach used in Refs [36l . 
[37j , we find the unitary transformation that diagonalizes 
the Jaynes-Cummings Hamiltonian H s 



D 



-A(JV,)/_ 



where 



A(iV 9 ) = - 



arctan (2\y/Nq) 



N„ 



a) a + II e 



(3.5) 

(3.6) 
(3.7) 



N q is an operator representing the total number of excita- 
tions, and n e is the projector on the qubit excited state. 
Applying this transformation to H s yields 



H 



D^D 

huj r a a H 



hu> a — ( 1 — 

2 2 



1 + 4X 2 N q I a z . 



(3.8) 



As it should, the eigenenergies of this Hamiltonian are 
the same as those presented in Ref. [l4[ if n is taken as 
the eigenvalues of N q and each eigenenergy is shifted by 
a constant Huj r . 

In this basis, the qubit is dressed by the field. As a 
result, qubit operators acquire photon part and similarly 
for field operators. For example, under the transforma- 
tion D, a z and a^a become 



2A 



y/1 + 4X 2 N q J y/1 + 4X 2 N q 
a z (XI+-a z /2) 



v/1 + 4A^A 2 



/+ (3.9) 



(3.10) 



Both these operators now involve the off-diagonal opera- 
tor 7j_. 



Expanding Eq. (|3 . 8[) to third order in A = g/A [one 
order up from Eq. (|3.3[) ] we find 



Hz 



where 



i h(u) r + C)a' a + h 
+ h({a^a) 2 a z , 



2 X [a ] a+- 



1 



1± 

2 



X = . 9 2 (l-A 2 )/A 



(3.11) 



(3.12) 



is the modified value of Lamb and Stark shift per photon. 
In addition to a correction to these values, the third order 
expansion yields a squeezing term (a'a) 2 of amplitude 
C = -.97 A 3 . 



C. The drive Hamiltonian under the exact 
transformation 

It is important not only to transform H s but also the 
drive Hamiltonian Hd- To do so, we first consider how the 
qubit and field ladder operators are transformed under D. 
Contrary to a z and N, the transformation does not lead 
to a compact result. To order C(A 3 ) for a and C(A 2 ) for 
cr_, wc find 



1 



AV- 



+ A 



1 - 3A 2 a\ 



1 



+ X A a z a+ 



(3.13) 



1-A 2 aW- 



Xaa z — A 2 a 2 <74 



(3.14) 



such that the drive Hamiltonian Eq. (|2.4p becomes 
X 2 a z . 



-iuJkt 



h.c. 



1 - 3A 2 f a f a+ - 



cr_|_e 



—iu) k t 



h.C. 



(3.15) 



With tok ~ w r , the first line of the above equation is 
responsible for measurement of the qubit. Due to the 
A term, the effective measurement drive strength is af- 
fected by the state of the qubit. It will be slightly larger 
or smaller depending on the qubit being in its excited or 
ground state. As will be shown later, this lead to small 
corrections to the ac-Stark shifted qubit transition fre- 
quency and measurement-induced dephasing rate. More- 
over, choosing u>k ~ w tt , one could take advantage of the 
second line of Eq. (|3.15[) to coherently control the qubit. 
Again due to a A 2 term, the effective strength of this con- 
trol drive will be modulated by the number of photons 
in the cavity. 
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IV. DISPERSIVE EFFECT ON THE MASTER 
EQUATION 

To obtain a complete description of the system in the 
dispersive regime, we also need to apply the dispersive 
transformation to the bath-system coupling. In princi- 
ple, this can be done by transforming the operators en- 
tering the dissipative terms of the Lindblad master equa- 
tion (|2.5[) . Once transformed, these terms will typically 
involve both qubit and field operators which correspond 
to probing the environment at different frequencies than 
the untransformed dissipative terms. Since the master 
equation is obtained in the Markov approximation, this 
frequency information is lost. 

Here, we go beyond this approximation by rederiving 
the qubit-resonator master equation. We first apply the 
dispersive transformation on the system-bath Hamilto- 
nian and then trace out the bath degrees of freedom to 
finally obtain a master equation in the dispersive frame. 



A. System-bath Hamiltonians 

Energy damping of the resonator (k) and of the qubit 
(7) can be modelled by coupling to baths of harmonic 
oscillators with free Hamiltonians [33| 



H 



B-f 



ub],(uj)b K (u))dw 
u)bt.(u))b~[(u>)du>, 



(4.1) 



where b\. (tjS) and & k . 7 (oj) respectively create and anni- 
hilate an excitation of frequency w in the resonator or 
qubit bath. Coupling to these baths is described by [33[ 



H K = ifi 



(/»&]>) - h.c.) (a + a?)ckv 



ih 



\J d 7 (w) (u>) - h.c.) <r x duj, 



(4.2) 

where di(u>) is the density of modes of bath i and fi(u>) 
represents the coupling strength of the mode of frequency 
cu to the resonator or the qubit. 

Dephasing in the bare basis occurs due to slow fluctu- 
ations of the qubit transition frequency. For example, in 
a superconducting charge qubit this is primarily caused 
by charge noise [H, HH ■ Dephasing can be modeled by 
adding the Hamiltonian 

H v = %vf v {t)a z . (4.3) 

In this expression, f v (t) is a random function of time 
with zero mean and v is characteristic of the magnitude 
of the coupling of the qubit to the fluctuations. Defin- 
ing f<p(t) = fp^e^du), H v can be written in fre- 
quency space as 



H v = hva z 



B. Dispersive master equation 

As shown in appendix[B] applying the dispersive trans- 
formation on the above system-bath Hamiltonians and 
integrating out the bath degrees of freedom leads to the 
master equation 

•D ■InD 1 ttD Dl 

Q = ~H H s + H d \ 

+ KV[a(l + A 2 ( 7 2 /2)]g D + lK V[°-](P 

+ [cr_{l - A 2 ( a ta + 1/2)}] g D + Kl V[a(j z }g D 

+ j v V[a z {l - 2A 2 (a f a + !/2)}}g D /2 



(4.5) 



where we have defined the rates k 
7 = 7 , k 7 = A 2 7 r with 



K rj 7k — A 2 /i a , 



/ v (w)e iut dw. 



(4.4) 



k p = 2ird K (uj p )\f K (uj p )\ (4.6a) 

7 P = 27tcZ 7 (cj p ) |/ 7 (w p )| 2 (4.6b) 

7 V = 2is 2 S(lu -> 0) (4.6c) 

7±A = 4AVS(±A). (4.6d) 

In obtaining these result, we have taken into account the 
fact that the spectral weight of the environment can be 
non- white. As a result, although we obtain a Marko- 
vian master equation, the rates depend explicitly on the 
qubit and resonator environments at different frequen- 
cies. As is explained in appendix [BJ in obtaining these 
results it was assumed that noise at the various relevant 
frequencies are independent. This assumption is valid if 
the noise is relatively weak and the various frequencies 
entering the expression for the rates are well separated 
one from another. For superconducting charge qubits 
which experimentally show long coherence times (up to 
2 /xs [10]) and in the dispersive regime (where A, and 
thus the frequency separation, is 10 GHz), the above 
model is accurate. 



V. EFFECTIVE QUBIT MASTER EQUATION: 
ELIMINIATION OF THE CAVITY 

In this section, we eliminate the resonator degree of 
freedom from Eq. (|4. 5[) to obtain a master equation for 
the reduced qubit density matrix in the dispersive frame. 
Building on Ref. [2l[, this is done by first moving to a 
rotating frame for the cavity, and then using a polaron- 
type transformation to displace the cavity field back to 
the vacuum. From this frame it is possible to consider 
only the two classical fields a e and a g . These fields cor- 
respond to the average value (a) of the cavity field if the 
qubit is in the excited or ground state. 

The resulting master equation is valid as long as the 
resonator state does not deviate too much from a super- 
position of coherent state. This can be formalized by two 
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requirements. The first is 

n <C n h 



K 
ICl 



(5.1) 



where n K is the ratio of the rate £ at which the non- 
linearity is squeezing the resonator state and the rate k 
at which these deviations are taken back to coherents 
state by damping. The second requirement is 



11, 11 < «, 



(5.2) 



where 7^ and 71 are given in Eqs (|5.7[) and (|5.8[) and are 
the rates at which the superposition of the coherent states 
a e and a g are getting mixed. This condition implies that 
the rate at which the superposition of the coherent states 
gets mixed is much slower than the rate of photon loss. 

As derived in appendix [C] in a frame rotating at ui m 
for the resonator, the effective qubit master equation is 



2 L j 2 
7 i P[a_] / , D + 7T I?[a + ]p D , 



(5.3) 



where p = Tr r (g) is the reduced density matrix of the 
qubit. The parameters introduced in this master equa- 
tion are 



lo™ =uj a ' + 2 [x + C(l + n e + n g )) Rc[a g a* e ] - C(n 2 g + n\) 

+ A 2 Re[e m/ i*] - 7 ' A+7 2 A " A27 Im[a 3 a;] (5.4) 
(n e + n g + 1) 



TVeff ~7y 



1 



2n cr ;t. 

T d =2[x + C(l + n e + n g )}lm[a g a* e } + X 2 lm[e m p' 
Rc[(7_ A a 9 - 7A"e + lX 2 a e )/3*] 



+ 



71 =7i 



7T =7-a™ 3 



2n K 



+ 7k +7A("-e + 1) 



(5.5) 
1 

(5.6) 

(5.7) 
(5.8) 



where the classical parts of the field a g and a e (pointer 
states) satisfy (i G {e,g}) 



A 2 



i < A' ± 



X + 2C n. 



k(1 ± A 2 ) + k 7 + 7±a - 5j^y 



(5.9) 



a,-, 



with the top sign for j = e and the bottom sign for j = g, 
and with AJ, m = A rm + £, where A rTO = aj r — cj m . In 
this expression, 5j. e is the Kronecker delta. Following the 
notation of Rcf . 2l| , we have used 



[3 = a e 



(5.10) 



Moreover, n e ^ g = \ot e ,g\ is the number of photons in 
the cavity when the qubit is in the ground or excited 
state. With A 2 = 0, the results of Ref [2l[ are correctly 
recovered. 

We now turn to a physically motivated description of 
these results. First, w Q D is the qubit transition frequency 
with Lo a ' — w a + X being the and Lamb shifted qubit 
frequency. The remaning terms are the ac-Stark shift. 
Then, Eq. (|5.5p is the qubit's pure dephasing rate. The 
first term is the bare pure dephasing rate, which now 
depends on the photon number n e + n g because of the 
dressing of the qubit by the field. Interestingly, this rate 
decreases with photon population as dressing increases 
with photon number and the photons are unaffected by 
qubit dephasing. This rate always remains positive as n e 
and rig must always be smaller than n cr j t . 

The term Td, which is defined in Eq. (|5.6p is 
measurement-induced dephasing. The first term comes 
from information about the qubit state contained in the 
frequency dependence of the pointer states [see Eq. (15. 9p ], 
the second term is the qubit information encoded into the 
driving part of the pointer states [see Eq. (|3.15p ]. and the 
last one is qubit information encoded into the decay rates 
of the pointer states. As in previous work 0, [2JJ, this 
decay rate can be negative. This is due to the recurrence 
of the qubit coherence which physically comes from the 
information about the qubit state which was lost into 
the cavity being transfercd back into the qubit. Positive 
constraints on the master equation bound how negative 
this rate can be but since our model is based on a phys- 
ical model which is positive in the enlarged cavity-qubit 
space these negative rates will never lead to an unphysi- 
cal state. 

Eq. (|5.7p represents the effective qubit decay rate. It's 
main contribution is proportional to 71 and, again, is 
reduced by dressing (but can never be negative). The 
second term j K ~ X 2 k is the Purcell effect which corre- 
sponds to qubit decay through the photon loss channel 
[34|. The last term of 7^ and 7-f are particularly inter- 
esting. They describe, respectively, additional relaxation 
and excitation of the qubit due to the photons populating 
the resonator. 

For the remaining of this section, we consider a purely 
white noise approximation, and a measurement drive 
at the resonator frequency, which ensures n = (n e + 
n g )/2 sa n e « n g . As a result of these contribu- 
tions, photons injected in the resonator for the measure- 
ment appear to the qubit as a heat bath of temperature 
T = (TiLOr/k-Q)/ log(l + 1/n). This effective temperature 
depends on the measurement drive amplitude and fre- 
quency. Under measurement, the qubit therefore suffers 
from additional mixing, something which can reduce the 
quantum non-demolition aspect of the readout. From 
Eq. (|4.6p . we can write l^/i = 7±a" ~ 2A 2 7 ¥ ,n and 
these mixing rates are therefore due to both photon pop- 
ulation and pure qubit dephasing. This can be under- 
stood in the following way. Let's assume the photon to 
initially be in the uniform superposition (jg) + |e))/\/2 
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FIG. 3: Dispersive energy diagram. The full line represent 
the bare states while the straight dashed lines represent the 
qubit-resonator dressed states. 



and the resonator in the vacuum state |0) (measurement 
drive is initially off). This initial state is schematically 
illustrated by the light gray dots in Fig. [H When the 
measurement drive is turned on, the photon population 
increases to reach a poisson distribution centered about 
an average value n. This is schematically illustrated by 
the black dots in Fig. [3] Because of the qubit-resonator 
coupling, the state |e,n) acquires a component \g,n + 1) 
and, likewise, \g,n) acquires a \e,n— 1) component (il- 
lustrated by the wiggly arrows). The amplitude of this 
qubit-resonator coherent 'mixing' increases with photon 
population as ~ X^/n. As illustrated in Fig. 7 of Ref. [Til ], 
in the absence of dcphasing this mixing is completely co- 
herently undone once the measurement drive is turned 
off leading to a QND measurement. However, in the 
presence of pure qubit dephasing, the phase coherence 
in the qubit-resonator dressed states can be lost, leading 
to effective downward and upward incoherent transitions 
between the qubit states. Given this, one should expect 
this rate to be proportional to the square of the qubit- 
resonator mixing amplitude X^/n and to the dcphasing 
rate, the result obtained in Eqs. (|5.7|) and ()5.8j) . 

In summary, the rates 7j and 7^ are due to dressing, 
by the resonator field, of the qubit operator a z causing 
dcphasing in the bare basis. We will therefore refer to 
this as dressed dcphasing. 



A. Numerical comparision with the full master 
equation 

To verify the validity of the previous results, we have 
done extensive numerical calculations in the limit n <C 
n K to compare results obtained from the reduced mas- 
ter equation (|5.3|) to those obtained from the qubit- 
resonator master equation (|2.5p . The results obtained 
from Eq. (|5.3p are also compared to those obtained from 
the linear approximation of Ref. [2l[. From this latter 
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FIG. 4: (Color online) Comparison between the exact master 
equation (|2.5[1 and the model (|5.3p . a) A typical time evo- 
lution of cr x (full black), G y (full grey) and a z for the exact 
result (full black), the non-linear (dashed blue) and the lin- 
ear (dotted red) models. The parameters and the initial state 
are the same as Fig. [2] b) Maximum of the trace distance 
for the linear (red squares) and the non-linear (green trian- 
gles) models, for j^/2ir = 0.5, 0.3, 0.1MHz (dotted, dashed, 
full lines). 



comparison, it will be apparent that the non-linear model 
obtained here is much more accurate, while adding essen- 
tially no additional complexity in numerical simulation. 

Figure QJi) presents a typical time evolution of the 
qubit as obtained by the numerical integration of the 
full master equation ()2.5p [full black line], the reduced 
model Eq. (|5.3|) [dashed blue line] and the linear model 
of Ref. [2l| [dotted red line]. The time evolution of (<r x ) 
and (a y ) obtained from these three models are indistin- 
guishable. However, because the linear model does not 
capture dressed dcphasing, only the non-linear model re- 
produces the correct equilibrium value of (a z ) 

We note that the numerical results obtained using the 
full master equation have been time-averaged to get rid 
of small amplitude fast oscillations. These oscillations 
are not contained in the effective models because of the 
various rotating-wave approximations that have been per- 
formed analytically. Experimentally, this averaging is ef- 
fectively performed due to the finite bandwidth of mea- 
surement apparatus. Moreover, for simplicity, for the 
numerical results, a white noise spectrum was assumed. 
We have therefore taken 7 K = X 2 k, k 7 = A 2 7i and 
7±A = 2A 2 7 ¥ , throughout this section. 

Figure |4p) shows the maximum of the trace distance 
(over time) as a function of measurement power for three 
values of the pure dcphasing rate 7 y . The red curves with 
square dots are the trace distances between the full and 
the linear reduced model, while the green curves with tri- 
angle dots are the trace distances between the full and the 
non-linear reduced model. Unsurprisingly, as the mea- 
surement power is increased, the trace distance between 
the reduced models and the exact solution increases. As 
shown the three different curves for both models, the dis- 
tance also increases as the dcphasing rate is increased. 
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However, the non-linear model obtained here is clearly 
much more accurate than the linear one, it captures the 
physics of dresscd-dephasing. The non-linear model also 
shows much less variation in the trace distance with de- 
phasing rate. It is worth pointing out that the maximum 
measurement power used in Figure SJa) corresponds to 
a very conservative photon population of the resonator 
n ~ 0.4n cr it, much lower than the critical photon number 
where non-linear effects were thought to become impor- 
tant BIU. 

The effective model developed in this section is both 
accurate and much less demanding numerically than the 
full numerical integration of the qubit-resonator Hamil- 
tonian. It should therefore be a useful tool to study the 
dispersive regime of circuit and cavity QED. 



VI. QUBIT POPULATION AND EFFECTIVE 
DAMPING RATE 



In this section, we focus on the dependance of the qubit 
mixing rate on photon population and dephasing rate, 
and on the steady-state value of (<t 2 ). These two quan- 
tities could be measured experimentally as a test of the 
present model. 



A. Photon number dependant qubit decay rate 

A remarkable feature of the non-linear model is that 
the qubit up, down and dephasing rates depend on the 
photon population. In particular, the effective qubit mix- 
ing rate is given by 



7eff(^ e 



,) = 7l0e S )+7T( n 9s) 



7i 



1 - 2X 2 



+ 7k 



(6.1) 



+ jA(n es + 1) + 7-A« S s, 



where rii S = \di S \ are understood as the steady-state so- 
lutions of Eq. (|5.9[) . Interestingly, in the situation where 
n es ?s n gs and for white noise, such that 7a = 7- a = 
2A 2 7 V , if 7 V < 7i/2, then increasing photon population 
leads to a decrease of the effective mixing rate. On the 
other hand, if 7 V > 7i/2, increasing photon population 
leads to an increase of the mixing rate. This is again a 
consequence of dressing of the qubit by the photon field. 



B. Measurement-induced heat bath 

From the reduced qubit master equation (|5 . 3[) . the 
steady-state value of (a z ) can be expressed as 



71 ~ 7T 
71 +7T 



= -1 



2 7tKs) 
7eff(rc es ,rv) 



(6.2) 



While the linear model would predict (er 2 ) s = — 1, the 
second order term 7_a causes a deviation of (cr z ) from 



this value which increases with n gs . This deviation is 
indicative of the breakdown of the QND aspect of the 
qubit measurement. 

When comparing expectation values, it is of course im- 
portant to compare expressions computed in the same 
basis. As a result, it is useful to transform (|6.2[) to the 
bare basis [46| . This is done by applying the dispersive 
transformation to er z , sec Eq. (|3.9[) . from which we obtain 



1 



'*/sb 



1 + 4A 2 



(oto), + 



(6.3) 



with (cr z ) s given by Eq. (|6.2|) . The last term of Eq. 
was neglected in the above expression as it oscillates 
rapidly in the rotating frame and the operator N q has 
been replaced by its average value. This expression for 
(cr z ) sh is consistant with Eq. (30) of Ref. [3| when notic- 
ing that there n corresponds to the total number of exci- 
tations while it is here the average number of photons. 

It is interestingly to note that, because of the asym- 
metry of the expression for (cr z ) with respect to n e and 
rig, the steady-state value of a z depends on the measure- 
ment frequency. Indeed, for a measurements of the phase 
where uj m = uj r , in the steady-state n es sa n gs sa n s and 
the effect will depend on the average number of photons 
in the resonator. On the other hand, for amplitude mea- 
surements with Lu m = Lu r + x, Ti es 3> Ti gs and the depar- 
ture from -1 should be less important. 



C. Comparison with exact numerics 

To compare the results of the analytical expressions 
Eq. (|6.3[) and (|6.1[) to numerical integration of the full 
resonator-qubit master equation (|2.5[) . we initialize the 
qubit in its excited state and the resonator in the cor- 
responding steady-state with a continuous measurement 
drive of amplitude e m and frequency tu r . In the absence of 
coherent driving at the qubit frequency, the qubit then 
simply decays to reach a steady-state value of a z . By 
fitting the time evolution of (a z ) as obtained from nu- 
merical integration of Eq. (|2.5|) to 



fit) = Ae^ + (a z ) 



sb ' 



(6.4) 



we extract the exact effective decay rate and steady-state 
mean value of <j z . 

These results are shown in Fig. the analytical ex- 
pressions (|6.3p and (|6.1[) (lines) in addition to the val- 
ues extracted from numerical solution of the full mas- 
ter equation (symbols) are plotted. With the parame- 
ters used here (see caption), the critical number of pho- 
tons is n cr j t = 400 such that the figure shows results for 
n/n clit . < 0.4. 

Figure [5^,) shows the steady-state value of (a z ) as a 
function of the measurement amplitude for increasing val- 
ues of the pure dephasing rate j v = 7±a/2. The bottom 
line (black stars) corresponds to "f v = 0, which in turns 
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FIG. 5: (Color online) a) Steady state value of (<r z ) and b) 
effective decay rate 7 e ff as a function of measurement power 
for various values of -y v /2iv = 0.0, 0.05, 0.2, 0.5 MHz (black 
stars, grey squares, red circles, blue triangles). Symbols are 
extracted from numerical solution of Eq. (|2.5[1 . Lines corre- 
spond to Eqs (|6.3p and (|6.ip . The numerical simulations were 
done with the same parameters as Fig. [2] with the qubit ini- 
tially in its excited state and the cavity in the corresponding 
steady-state with a continuous measurement drive of ampli- 
tude e m . The top x-axis is the approximate average photons 
number in the cavity. 



corresponds to the effective heat bath being at zero tem- 
perature. As a result, in the dispersive basis, (cr z ) s = — 1 
and the deviation from —1 is only caused by the change 
from dispersive to bare basis. The lines lying above this 
result correspond to 7^/2^ = 0.05, 0.2, 0.5 MHz (dashed 
grey, dotted red, dashed-dottcd blue). Clearly, even for 
a relatively low number of photons compared to n cr ; t ., 
mixing of the qubit excited and ground states by the ef- 
fective heat bath can be significant if the pure dephasing 
rate is large. 

While the reduced model is extremely accurate for 
"f v = 0, it always slightly over-estimates (cr z ) b for non- 
zero dephasing rates. Since the error is always positive, 
this can be interpreted as being due to the effect of a 
higher order terms in the dispersive approximation. Nev- 
ertheless, the analytical model is at most ~ 10% away 
from the exact numerical results for the range of param- 
eters shown in Fig. [5] 

Figure O)) shows the effective decay rate 7 e s as a func- 
tion of measurement power for the same dephasing rates 
7 V as in pannel a). The lines correspond to Eq. (|6.ip 
while symbols are extracted numerically. As expected 
from the discussion surrounding Eq. (|6.ip , if the dephas- 
ing rate is negligible (black stars), the qubit effective de- 
cay rate falls below the bare decay rate 7i/27r = 0.1 MHz 
as the measurement power is increased. However, for 
2j v > 7i (red circles and blue triangles) the effective de- 



cay rate increases, again as expected from the model. In 
this latter case, the photon number and dephasing depen- 
dent qubit mixing simply overwhelms the aforementioned 
decrease in j e g which is no longer visible. For interme- 
diate dephasing, 27^ = 71 (gray squares), these two pro- 
cesses cancel each other and the effective decay rate is 
almost independent of the measurement amplitude. 



VII. DISPERSIVE EFFECTS ON THE 
QUANTUM TRAJECTORY EQUATION 

The master equation description of the dynamics does 
not take into account the result of the measurement. To 
include this information, we use quantum trajectory the- 
ory [H, H3, IH, EH and derive the evolution equation 
for the conditional state, or quantum trajectory equation 
(QTE). This was already study in Ref. [2l| for the linear 
dispersive model and is extended here to incorporate the 
non-linear effects. 

When monitoring the resonator bath, characterized by 
the rate k, an observer would in principle see two sig- 
nals: one at the qubit frequency and one at the cavity 
frequency (47[. To derive the QTE, it is assumed as above 
that the relevant bath frequencies for the resonator bath 
are well separated such that they can be treated as two 
separate Markovian baths with relevant frequencies u> r 
and u) a . That is, for an infinitesimal interval dt the uni- 
tary operator describing the resonator bath is given by 
Eq. Ip9a]) . 

In a homodyne measurement, with a local oscillator 
set to ~ u r , the bath is projected in an eigenstate of 
the operator dB K , r + dB Kr (where dB K ^ r is defined in 
appendix [B| . with measurement result J [33]. For many 
such measurements, each separated by a time dt and with 
result Jk, the state conditioned on the complete record 
3(t) = { Ji, Jk} can be expressed as HU 

Qj(t) = ej(t)/Pr(J), (7.1) 

where Pr(J) is the probability for observing record J 

Pr(J) = Tr[gj(t)] (7.2) 

and gj (t) is an unnormalizcd conditional state given by 

= E Mj k , Fk (dt)...M Jl>Fl (dt)g(0)M\ i Fi (dt) 

F 

(7.3) 

Following the notation of Ref. (43|, Mj fc> p k is the Kraus 
operator of the k th measurement and Fk = {/2k, fsk, fik] 
the results of fictitious measurements performed on the 
baths coupled to the operators L2, L3 and L4 defined in 
Eq. (jBip . Since these latter measurements are fictitious 
and the conditional state does not depend on them (we 
sum over all fictitious results to obtain gj), it is possible 
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to choose the fictitious observable at will [43j. For sim- 
plicity, we assume fictitious homodyne measurement of 
the unobserved operators Lj with j = 2,3,4. Using the 
evolution operator Eq. (|B8[) . this leads to the following 
Kraus operator 

M JtF (dt) = (F, J\ U(t + dt, t) |0, 0) 

= V /t J,f[1 - iH tot dt + ^fnLxJdt + K,L\L\dt/2 
+ F T Ldt - VLdt/2] 

(7.4) 

where L = { v /k^X 2 , y/j^L 3 , % / 7 jvL 4 } and Tj,f is the gaus- 
sian probability measure 



appendix iDl this reduced QTE takes the form 



Tj tF dJdf = 



{2ir/dt)'- 



exp[— ( J 2 + F 2 )dt/2]dJdF. (7.5) 



For continuous monitoring, the time step dt between 
measurements tend towards 0. In this limit, Eq. (|7.3p 
leads to the QTE whose ensemble average is the uncon- 
ditional master equation. The corresponding QTE, in 
Ito form, for the measurement operator Eq. (|7.4p is 



=£ D ^ + 2^M[I 4> (1 + \ 2 a z /2)] 6 ^m 
+ i^j[Q 4> (l + X 2 <J z /2),g^(t) 



(7.6) 



where CP is given by Eq. (|4.5[) . In this expression, we 
have defined the ^-dependent field quadratures 21^ = 
ae-^ + aV* and 2Q = -iae'^ + ia^e i(/> . Moreover, 
the superoperator M. [c] is defined as 

M[c]q =(c- (c) t )g/2 + g(c - (c) t )/2, (7.7) 

where (c) t = Tr (cg>y(i)) and the measurement outcome, 
J can be expressed as 



J(t) = 2^j(l^l + \ 2 a z /2)) t +f(t), 



(7- 



where £(t) is Gaussian white noise and r\ is a detection 
efficiency parameter included for completeness (2~i| . 

While only the signal at the cavity frequency was taken 
into account here, it is interesting to point out that the 
signal at the qubit frequency could also be measured. 
This was done experimentally in Ref. [26| to perform 
qubit state tomography. However, the signal at that fre- 
quency is in general much weaker than the signal at the 
cavity frequency. As a result, while by itself the former 
signal would lead to a very inefficient QTE (which would 
have a similar to that of a direct homodyne measure- 
ment of the qubit) this additional information could be 
included in the present treatment to realize even more 
efficient qubit measurements. 



A. Effective qubit quantum trajectory equation 



Using the polaron transformation Eq. (|C1[) . it is possi- 
ble to obtain a reduced QTE for the qubit. As show in 



p° = CP p v + ^~if)M[<j z ]p^(t)(J(t) ^/f~(t){a z ) t ) 
_.^/KJf)_ ^ j/0 D (t)] (J(t) _ y /f^)(a z ) t ), 

(7-9) 

where r c ;(£) is the rate at which information comes out of 
the resonator and T^^it) represents extra non-Hciscnbcrg 
back-action from the measurement. J(t) is the processed 
record coming from the resonator and is given by pll | 

j(t) = ^/f~(a z ) t + £(t). (7.10) 
This quantity is linked to the homodyne current by 



J(t) = J(t) + ^rj \n(t)\ cos((9 M - 0) 
_A 2 

rc^y I^WI cos (°i3 - 



(7.11) 



Eq. {779]) has the same form as the QTE found in Ref. [21| 
for the linear model, appart from second order correc- 
tions to the rates T c i(t) = r]T m cos 2 (9 m ) and Tba(i) = 
rjF m sin 2 (# m ), where 



1 |^|cos(6>3 - 6> M ) ImI 2 

4 |/3|n C rit. 64|/?| 2 n c 2 r 



4> - 9 /3 + Im < In 



8 1/31 n crit . 



(7.12a) 
(7.12b) 



These expressions are valid for <fi — 613 G [0, 7r/2], with 
Op = arg(/3), = arg(/i) and (3, /i defined in Eq. (|5. 10[) . 
The second order corrections do not change the physics 
in an important way since, as in the linear model, it is 
possible to choose the phase (f> of the LO optimally such 
that Tba(i) is zero. The corrections have the effect of 
reducing r c j(£) in comparison to what is obtained in the 
linear model (2l| . 

To demonstrate the different features of the reduced 
QTE, Fig. [S] presents a typical trajectory for 3 differ- 
ent measurement powers using the same parameters as 
in Fig. [21 As in the linear case, increasing the measure- 
ment power localizes the qubit state on one of its basis 
states. As a result, although the QTE is based on homo- 
dyne measurement, we do not expect diffusive but rather 
jump-like trajectories. Moreover, because of the effec- 
tive upward rate 7^ which increases with measurement 
power, the trajectories show telegraph noise rather than 
a single jump to the ground state. These predictions can 
be experimentally tested once single shot measurement is 
achievable. From these results, the waiting time between 
jumps can be compared to 77 and y±. 

The SNR can be defined as 



SNR 



Td 

7eff 



(7.13) 
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FIG. 6: Typical trajectories for 7^/2^ = 0.5MHz and mea- 
surement amplitudes of e m /27r = MHz (bottom), e m /2n — 
10 MHz (center) and e m /2?r = 20 MHz (top). The other 
parameters are the same as in Fig. [2] The initial state has 
(<r z ) = 1 with zero photons and a measurement drive starting 
shortly after t = 0. 



have obtained an effective stochastic master equation for 
the qubit. In the single-shot limit, this equation predicts 
that the measurement-induced heat bath should lead to 
telegraph-like jumps in the measurement response. More- 
over, the non-linearity have been shown to lead to a reduc- 
tion of the expected signal-to-noise ratio, a result quali- 
tatively consistent with experimental observations [44| . 

There are various ways to test experimentally these pre- 
dictions. First, the qubit effective decay rate 7j (n)+j\, {n) 
and the steady-state value of (a z ) can be measured ex- 
perimentally and compared to the results obtained here. 
Second, when single shot measurements become possi- 
ble in circuit QED, the effect of the upward transitions 
caused by the measurement induced heat bath should be 
observed. The waiting times between the upward and 
downward transitions can then be related to the rates 
7l(n) and 7j,(«) obtained here. 



which at the optimal point is 7?r m /7 c ff. As in the linear 
model, r m is proportional, through \f3\ , to the cavity 
pull x + C a<!a times the number of photons in the cavity. 
However, and contrary to the linear model, the cavity 
pull decreases with increasing photon population. More- 
over, we have shown that the measurement photons act 
as a heat bath with j c ff oc n. Therefore, unlike the linear 
predictions, which predicts the SNR to increase linearly 
with photon population, we expect the SNR to saturate 
at higher photon numbers. This is discussed in more de- 
tail in Ref. [44|. It is important to point out that the 
main contribution to this effect is the reduction of the 
cavity pull, and not the measurement-enhanced mixing 
rate. 



VIII. CONCLUSION 
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APPENDIX A: EXACT DIAGONALISATION OF 
THE JAYNES-CUMMINGS HAMILTONIAN BY 
UNITARY TRANSFORMATION 



We have investigated circuit QED in the dispersive 
regime. To take into account large photon population 
of the resonator, useful for qubit readout, we have shown 
that it is necessary to push the dispersive treatment to a 
higher order. We have done this while taking into account 
the effect of dissipation and external microwave driving. 
In particular, we have obtained a Markovian model for 
the effect of dissipation that takes into account frequency 
dependence of the environment. 

Building on our previous work [2lj . we have then 
traced over the resonator states to obtain an effective 
master equation for the qubit valid in the limit n <C 
{n K = k/C, n cr i t = A 2 /4g 2 }. A striking feature of the re- 
sulting master equation is that the qubit relaxation and 
dephasing rates now depend on the number of photons 
populating the resonator. Moreover, in the presence of 
pure dephasing, we have shown that measurement will 
cause excitation of the qubit. In other words, the photon 
population of the resonator act as an effective heat bath 
on the qubit. This can lower the effective 'quality' of the 
dispersive QND readout of the qubit. 

Finally, using the quantum trajectory approach, we 



In this appendix, the Jaynes-Cummings Hamiltonian 
is diagonalized exactly using a unitary transformation. 
In order to simplify the notation, we introduce the com- 
mutation linear application 



C A B = [A, B] 



m times 

C%B=[A, [A,[A. 



.,B}]]. (Al) 



In term of this superoperator, Hausdorff's relation can 
be written as 



z A Be~ A 



71 = 



-{IB- 



(A2) 



In the same way as for the linear approximation discussed 
in section llll Al the antihermitian operator /_ will be key 
in this diagonalization. Since it commutes with both H s 
and J_, another important operator is the total number 
of quanta N q defined in Eq. (|3.7|) . Moreover, with a uni- 
tary operator of the form 



D 



,-A(iV 8 )I_ 



(A3) 



where A is a function to be defined, N q can be considered 
as a scalar when applied on H s . 
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Before transforming H s using D, it is useful to intro- 
duce some important commutators. First, it is simple to 
show that 

C I _H a = hAI + . (A4) 
Using this result, transformation of H s by D yields 



with 



CaLA- = (-4)"A 2lI iV™/ + 



i2n+l 
"AZ_ 



(A5) 

(A6) 
(A7) 



Using this last result, we find 



H? = H n + h 



2hN q a z 




g sin (2Ay/TQ A (l - cos (2A-/AQ) 



4iV„ 



To complete the diagonalization, we take 



HN q ) 



— arctan (2\^JNq 
2JK 



(A8) 



(A9) 



such as to eliminate the off-diagonal term proportional to 
I + . Using this result, we finally obtain the exact diagonal 
form 



TiA 



H? = H - — 1 - Jl + 4X 2 N g a 



(A10) 



Using this result, we define the Lamb and ac-Stark 
shift operators as [we use HP(cva,a z )] 



S L = H?(0, 1) - H? (0, -1) - hLO a 
hA 
~ 2~ 

S s (a f a) = H? (a f a, 1) - H? (a f o, -1) - S L - huj a 
_ HA 



(l - Vl + 4A 2 ) 
Ja,l) - iJ s D (a f a 

(J 1 + 4A 2 (at a + 1) + Vl + 4A 2 ata 



(All) 



-1 - v / l + 4A 2 j . (A12) 

Developing these expressions in powers of A, wc obtain 

6 L ^h x + 0(\ 5 ) (A13) 
6 s ((Ja) « T^a + TiCC^a) 2 + C(A 5 ), (A14) 

with x = <A(1 — A 2 ) and £ = — g 4 /A 3 . These approxi- 
mate results are used in Eqs. (|3.3[) and ()3.11D . 



APPENDIX B: OBTAINING THE DISPERSIVE 
MASTER EQUATION 

1. Qubit relaxation and photon decay 

In this appendix, we find the effect of the dispersive 
transformation on the non-unitary part of the master 
equation. Applying the dispersive transformation on 
the Hamiltonians Eqs. (|4.2[) . moving to the interaction 
frame defined by the transformation exp[— i(H s + Hbk + 
HB-y)t/h] and performing a rotating-wave approximation 
(RWA) yields 



H? = ih [Li4(*,w P ) + La4(*.WB)] +h.c. 
H? = ili [L 3 z\(t,u a ) + L 4 z\(t,u r )] +h.c, 



(Bl) 



withLi = a(l + A 2 ^/2), L 2 = Aer_, L 3 = <r_[l-A 2 (a t a+ 
1/2)] and L4 = Xaa z . The bath operators 2j(t, u> p ) are 
given by 



Zi(t, U)p) 



uj p — Bi !P 



(B2) 

We have kept in Eq. (|B1|) only terms that will contribute 
up to order A 2 in the master equation. In obtaining this 
expression, we have taken the (dispersive) system Hamil- 
tonian as H s r* hu> r a^a+?iui a (T z /2, where io a should be un- 
derstood as the Lamb and ac-Stark shifted qubit transi- 
tion frequency and uj r should be understood as the cavity 
frequency shifted by the non-linearity (i.e. cu r + (). More- 
over, to perform the RWA, we have made the standard 
assumption that the system-bath interaction is limited to 
a small band of frequency Bi jP around the frequency tu p 
of the corresponding system operator Li, Bi tP <C ui v [331 ] . 

We now assume that, within the bandwidths J5i iP , the 
coupling constants fi(u> p ) and the density of modes di(uj) 
do not vary significantly. In this situation, the above 
system-bath Hamiltonians can be rewritten as 



H? = ih^Lxbl^t) + ih^jr a L 2 b{ a {t) + h.c. 
H? = ihy/^Lzbl, „(t) + ihyffrUbl ' (t) + h.c. 



(B3) 



where the decay rates are given by Eq. (|4.6[) and the bath 
temporal modes are defined as 



k lP (t) 



1 



27T J^p-Bi, 



dwbi(u>)e 



(B4) 



Since bi(u), bj(u>') = Sijd(u> — u>'), the commutator 

of two temporal modes is (after a change of integration 
variable) 



[M*),&L(0] 

hi 



2tt 



dujduj' ' 8{lo — u> — lo p + uj q )e 



-iut „ — icu't' 



(B5) 
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If we now take \u p — u> q \ 3> -Bj, p , -Bi,g for p =^ q then the 
above becomes 

K P (t)ib\ q {t>)] = i^e-Mt-O. ( B6 ) 

In other words, we assume the bath operators to be in- 
dependent. In the dispersive regime, this is a reason- 
able assumption since |oj p — | ~ A, where the detuning 
|A| > g is large. 

We finally make the standard and reasonable assump- 
tion that dissipation is not too strong, such that the time 
scales set by the decay rates k p and 7 P are much longer 
than the cutoff time 1/B^ p . In this situation we can 
effectively take the limit B^ p — > oo. This corresponds 
to the standard Markov approximation (33[, which was 
already successfully applied to describe circuit QED ex- 
periments [H EE El IM EE Ezt In this situation, the 
above commutation relation reduces to 

[kAt),b{ q (t')] = S itj 6 Ptq 6{t-t'). (B7) 

In this Markov, or white noise, approximation the evo- 
lution operators can be written in Ito form as 

U(t + dt,t) = U K (t + dt,t)U y (t + dt,t)e- iHt ° tdt (B8) 

with 

U K (t + dt, t) = cxp {-iy/K^LrfB^ - L\dB K , r ] 

-i^JT a [L 2 dBi a - L\dB K , a \}u K (t) 

(B9a) 

Uj(t + dt,t) = cxp <^-i^/j^[L 3 dB^ a - L\dB 1 ^ a ] 

-i^.\[UdB^ r - L\dB 1<r }} U 7 (t), 

(B9b) 

where dBi iP = b^ p dt is a quantum Wiener increment [33| . 

We now take the bath to be in the vacuum state and 
uncorrelated to the system at time t = 0. By tracing over 
the bath and keeping terms of order 0(dt) using Ito cal- 
culus, we obtain a Lindblad form master equation for 
the resonator-qubit system. In this master equation, the 
photon bath k leads to the damping superoperators [33[ 

KD[a{l + A 2 a 2 /2)]£ D + 1k V^-]0° (BIO) 

while the qubit bath 7 leads to 

7P [er_{l - A 2 (a f a + 1/2)}] g D + K 7 V[aa z ]g n . (Bll) 

These terms are the second and third lines of Eq. (|4. 5|) . 

2. Qubit dephasing 

For dephasing, we start with the Hamiltonian (|4.4|) . 
Moving to the dispersive basis, it becomes 

/oo 
U{uj)e lut doj, 
-OO 

(B12) 



where we have used the second order expansion of 
Eq. (|3.9[) . Moving to a frame rotating at the qubit and 
resonator frequencies we find 

Hf cp = hvo z {\ - 2X 2 N q )f (t) - 2hv\ct<r-Mt) 
— 2frh'\a<j+f-A{t), 

where 

/oo 
Ui^-^dco. (B14) 
-00 

The main contribution to dephasing comes from a small 
frequency band Bq centered around the frequency luq . In 
this situation, the integration boundaries in / Wo (i) can 
be reduced to 

fvott) = / Uicoy^-^du. (B15) 

For this rotating-wave approximation to be valid, it is 
required that B <C uiq [33[ . 

The Wiener-Khinchin theorem can be used to relate 
} v {lo) to its noise spectrum S(u>) [H[ 

E[U(uj)U(-uj')} = 8{lo - u/)S(u>), (B16) 

where E[-\ is an ensemble average. This allows us to write 
the uj component of the noise as 

= v^Re^), (B17) 
with £(u>) white noise obeying E[£(uj)] = and 

£K(«)M = <j( u -«'). 

Using these results, we now make similar assumptions 
as in the last section and take the noise spectrum S(u) to 
be constant within the small frequency band Bq around 
u>o. After a change of integration variable, / wo (i) can be 
written as 

f U0 (t) = / ° £(w + u )e iut du. (B18) 

J -B 

In this Markov approximation, we will again assume the 
noise spectrum to be relatively weak which implies that 
the time scale corresponding to dissipation is much slower 
than l/So [33] [. In this situation, we take Bq — > 00 which 
allows us to write 

f U0 (t) = v^o)U(*) ( B1 9) 

such that the transformed dephasing Hamiltonian be- 
comes 

Hjj> ep = fwy/S(jSja,(l - 2\ 2 N q )£ (t) 

- 2hvy / S(A)Xa^a^ A (t) (B20) 

- 2Tivy/S(- A) Aaa + ^_A (t) . 

The three £ wo (t) white noise terms in the above ex- 
pression now correspond to independent noises, centered 
around three different frequencies. 
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The above Hamiltonian leads to the following supcrop- 
erators in the resonator-qubit master equation 

l v V[a z {l - 2A 2 (a t a + l/2)}]g D /2 



(B21) 



with the rates given by Eq. (|4.6[) . These terms correspond 
to the fourth and fifth lines of Eq. ()4.5|) . 



APPENDIX C: THE POLARON 
TRANSFORMATION 

Following the approach developed in Ref [2l|, a re- 
duced master equation for the qubit is obtained in this 
appendix. To do so, we start from the dispersive master 
equation (|4. 5[) and go to the rotating frame defined by 
R = cxp[ia; rn a^at]. We then go to a frame defined by the 
polaron-type transformation 



U e D{a e ) + n g D(a g 



(CI) 



where D(a) is the displacement operator and a g ( e ) satisfy 
Eq. (|5.9[) . In the polaron frame, the field a is described 
by a classical part given by the complex variables a g and 
a e , and a small quantum part corresponding to quantum 
noise. 

The action of P on various system operators is given 

by 



pt a p 


= a + U a 






(C2a) 


PVaP 


= era + a'lla + aTI* 


+ |n Q | 2 


(C2b) 


Ptcr P 


= a_D^(a g )D{a e ) 




(C2c) 


pV,p 


= (?z 






(C2d) 


P t (a + a) 2 P 


-|n Q | 4 + 


"(2|H Q | 2 - 


-l)n*a + h.c. 






+ |n Q | 2 (4a 


f a+l) + 


(aaU* a 2 + h.c.) 


(C2e) 




+ (2a t aan* + h.c.) 4 


- (at a) 2 





where we have defined the projection operator 
U a = a g U g + a e n e , 



(C3) 



with A' rm = A rm + (. Taking into account the time- 
dependence of P, the transformed Hamiltonian reads 

H r>p = h dp _ ( in fi aa i + h c -) + ai m [n Q n*], (C5) 

the bar on the superscripts indicating that time depen- 
dence of the transformation is taken into account cxplic- 
itly. 

We also apply this transformation to the dissipative 
terms of the dispersive master equation (|4.5|) . For the 
first term of second line (k term), keeping up to order A 2 , 
we get 



a p I 1 



p = V 



a 1 



A 2 (7 2 



+ [<j z ,p] — (/3 + An a ) + h.c. 
+ -M\ 2 + X 2 (n e -n g )]V[a z ]p 



.Im[a 9 a*] 

1 2 
.1 



(C6) 



i- [-ill a (1 + X 2 a z ) at + h.c, p] + O (A 4 ) . 



V 



We then get for the 7 and k 7 terms 
J a- 



1 - 2A 2 1 n P 



1 

2, 



- a e \ 2 D(f3)a- (a f p + pa^a+D^ (J3) + h.c. 
+ 2a e A 2 pa^<j + a- + h.c. 

- A 2 {D^a-pa+a^ aD\/3) - pa'aa+a-} + h.c. 
-iX 2 [ia e n e at +h.c.,p] +C(A 4 ), (C7a) 



n 3 + \a e \ n n e . Using these results, we V[a p a z }p = V\aa z \p + ^-V[<r z ]p 



with |II a | n = |( _ 

apply the transformation P to the Hamiltonian H u = 
H? + H? to obtain 



4 2 2 



X 2 a 2 



= hA' rm \n a \ 2 + K\u a t 

+ n(n* e m + h.c.) (l + 

+ ft[^a + x + 2(x + C)|n Q | 
+ h A' rm n Q + ( x + OT\ a a, 

+ e m (l 

+ fi ' 



(C4) 



\ajl a pa) a z — patjM + h.c. 
i- [-io t n Q +h.c.,p] , 



(C7b) 



2 

2c|n c 

f h.c. 

K m +(x + Ca f a + AC\n a \ 2 )a, 



and, for the ~f v term, 



V 



1 - 2A 2 ( a) a + 



a) a 



?iCa t a t n 2 + h.c. + 2TiC,a) (J aX\ a o z + h.c. 



[1-2A 2 (n e + n 9 + 1)] V\a z \p 

- 2A 2 [o t aX'[(7 2: ]p + h.c] 

- 2A 2 [(V[a z ]p)a^n a + h.c] 

- 2A 2 [a^U a V[a z ]p + h.c] + O (A 4 ) 



(C8) 
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Finally, for 7 ±a, we have 



1. Reduced master equation 



T>[aro%]p = n g V[a%]p + V[a%a]p 

+ a g (D^(f3)a + pa^a t D(p) - patllg) + h.c. 

-^[-ia.atn.+h.c.p] (cga) 

V[a^ P <j p ]p = n e V[a p ]p + V[ata)]p 

+ a e [D{P)a)(j-po + D\f3) - pa)li e ] + h.c. 

- i^- [-ia e a^U e + h.c., p] . (C9b) 



The last line of Eqs ((C6J), §<JB) and ((C9J) act like a drive 
Hamiltonian. Wc will be able to cancel them with a e and 
otg given by Eq (|5.9[) . In the above expressions, the quan- 
tities n g = \a g \ 2 and n e = \a e \ 2 are the number of pho- 
tons when the qubit in the ground or excited state, and 
we have 2?[ct?] = V[D((3)<rJ\ and V[a^} = V[D^{j3)a+]. 

If we put all the results of this section together, we can 
write the polaron-frame master equation, which is given 
by applying the polaron transform on Eq. (|4.5[) . The re- 
sult is given by combining the results from Eqs (|C5IIC9[) 



•DP 



ff Dp ,g DP 



DP 



1-X 2 



t 1 
aa+ 2 



DP 



K 7 V[a p o z \g 



DP 



-V 



1 2A 2 f oW 



DP 



■ 7A2?[a t <J P ]g DP + 7 -A^[a P ^]? DP 



(CIO) 



In this section, we trace the transformed master equa- 
tion Eq. fC10|) over the resonator states to obtain an 
effective master equation for the qubit only. This is done 
by first expressing the total density matrix in the polaron- 
transformcd frame as 



DP 



DP 



^ ] ^ ] Qn,rr 
n,m=0 StS '£{e,g} 



\n,s)(m,s'\. (Cll) 



Since our goal is to obtain the effective equation in the 
original non-polaron transformed frame, we write the re- 
duced qubit density matrix in this frame as 



p D =Tr r (P, DP pt)= £ 

s,s'e{g,e} 



(C12) 



with 



Ps,s — 00,0, 31 Pe.g ~ ^n,m,m,n' (C13) 

where we have defined 



n,m— 



DP rr {pi DP 

Qi,j,s = Tr r [a 1 a g SiS 



(C14) 



with {s,s>} G {g,e}, \% P tM = Q™, e , g d p , q e-^<\ 
and d Ptq = (p \D[{3] \ q) is the matrix element of the dis- 
placement operator in the number basis. 

To obtain the master equation for p D , we simply find 
the equation of motion for the matrix elements of £? DP . 
More precisely, we will look at the equation of motion for 
P DP 



•DP 



,DP 

n . in a 



■ 2i((n - m)p° P lim+M 



-i [{A' rm + X + 4Cn e )(n - m) + ({n 2 - m 2 )] g: 

- iC«e [ 2n Qn-l,m+l,e + n ( n ~ l )Qn-2,,n,e] + W+l,m-l,e + m ( m ~ l )QnJn-2,e] 

- 2<"e [(2n - m)g™ P l+l e + n(n - l)e°-l,m,e] + 2i C a *e [i 2m ~ n )Qn+l,m,e + m ( m ~ i^m-l.e] 

- |[k 7 + k(1 + A 2 ) + 7a-2 7 A 2 ] 



n + m 



1 - 2A 2 in 



DP 

+ Ik + jAn e > g n „, e 



(C15) 



V 2 

■7-A 



/7A „ ,2\ i DP i * DP \ / n \2\ DP 

{— - 2 1 X ) (u e g n m+l e + a e g n+1<m>e ) - ( 7A - 2 7 A )g n+hrn+1 . e 
Tr r (D{P)a? m a n D\p){n g g™ + a g gf p a^ + a* g agfj + agf p a^)) 



From this equation, we see that the only way the clement 
Qo > fi,e depends on the other elements is through the two 
last lines. Moreover, the only way the elements n, m ^ 
can be populated from an element i < n,j < m is 



through the second, third and last lines. The rates at 
which these mechanisms act are of the order £ |a e | 2 = Qn e 
and 7 | = 7-A^g- On the other side, these elements decay 
more quickly than the 0, element because of the k term, 
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which we assume is dominant compared to k 7 , 7a and 
7A 2 . If the conditions n e -C n K = k/C, and 7j -C k are 
satisfied, we can assume there is no significant population 
of the n, m ^ matrix elements. We will have a similar 
equation for „, with the conditions being n g -C n re 
and 71 < k, where 7^ is defined at Eq. (|5.7p . If these 
conditions arc fulfilled, we can reduce the above equation 
and that for the g component to 



■D 

Pg-g 



D 

■liPe,e- 



(C16a) 
(C16b) 



On the other hand, the off-diagonal elements of the re- 
duced qubit density matrix involve off-diagonal elements 
of the resonator density matrix and we must consider the 
equation of motion for all the terms X^^n p q 

X DP - n DP rl o-iHorf-iflTmU rv*l X DP 

A n,m,p,q ~~ U n,m,e,g w P,q t ' lu t ±i ^l Ll g Ll e\ /y n,m,p,q 



+ PVP^n,m,p-l,q P V^l\i,m,p,q-1 (C17) 



2 " <■ v~ ' ) n,m,p,q' 



where the bar means that the state is not normalized and 
the linear measurement superoperator is 

M[c\g D = {cg D + g D c)/2. (D2) 



Moving to the frame defined by Eq. (|C1[) yields 

X 2 a 



-DP _ f DP -DP 

8j — L - Qj 



' KT] 



a 1 



e^QT 



I KT] 



where we have defined 



Re[/3 ]A^]^T +Re[A ]^ 

— K, Qj \ J, 



DP 



J 



(D3) 



A 2 /T 
fI+ — 1 



(D4) 

As it should, for A 2 = the three equations above are of 
the same form as those obtained in Ref. [2l[ . 

As before we now find the equations of motion for the 



Coefficients Q^e, Qn,rn,g> and A n,m,p.«r For thc 

clement, we find 



DP 

n,m,e 



If we do this and compute Qnme,g according to (|C10j) . 
we get an equation that can be reduced only to the cle- 
ment Qqq e g in the conditions stated above (n e , n g -C n K ). 
Considering that only the 0, element is ever populated 
significantly, the equation of motion is then 



Qn F m , e = (EES) + V^(Re[A ] + Re[/y )e°£, e J 

-i<t> rDP 



y/wj ( 1 + ~ ) [ e ^Qn+l.m.e + e *^f?n,fn+l,e] ^ 1 

(D5) 



■ DP 



Pe,g — A 0, 0,0,0 



7t 



= -i&a + djm[a g a*])p° g 
7J. ( x> 1 



7v e ff 



with 



w a + X + Re 



3 t (/?/?* 



AV 



D 



(C18) 



7y e ff. 



- C(n e + n 2 ) + (k + K 7 )Im[a g a* 
: 7 V [1 - 2A 2 (n e +n g + 1)] 

+ (k+K7 )^! + " a2 k-^) 



7a 
2 ' 



(C19a) 



(C19b) 



Using the expression (|5.9p for o; g ( e ) , we can combine 
the equations of motion for the reduced qubit density 
matrix p° to find the reduced qubit master equation (|5.3|) 
with the frequency and rates given by Eqs. (|5.4j) . (|5.5j) . 
d577l l and (15^1) . 



APPENDIX D: THE EFFECTIVE QUBIT 
QUANTUM TRAJECTORY EQUATION 

The QTE is derived using linear quantum measure- 
ment theory [HI, SI ■ The linear form of Eq. ([7T6]) is 



7^MlU(l + \ 2 a z /2)}Q°J 



+ i^T 1 [Q 4> {l + \ 2 a z /2),Q 1 ?] J, 



(Dl) 



with a similar equation for the g^^n „ component. For 



the A„ p g component, we find 



:dp 



(Rc[// ] +iIm[/3 ])A 



DP 

n,m,p,q 



A 



:dp 



1 + y ) \/n + le i0 A* +] 



1 - y ) v / ^TT e ^A^ +1 



J. 



(D6) 



In these expressions, which are the contribution of the 
Linblad term £ DP gij )P of Eq. (|E>3|) , the equation numbers 
refer to the RHS of the corresponding expressions. 

The added measurement and back-action operators in 
the evolution equations does not change the approxima- 
tion used in the previous section. Therefore, in the same 
limits, we can consider that the only relevant components 
arc ^cPcTe) ^o > o > o an d ^ao^o.o- We can then write 



Pe,e 



-D 

Pg,g 



Kg 



',3 

(ICTHaTl 



A ; ' 



Re[/} ] + Re[/3 ] J (D7a 



(ICT6bl) 



2 

K 7 / 



A 2 

^(1-y 



(ICT81) 



Ro[/i^] - Rc[/3^]J J (D7b) 

-i0-DP , # -DP ] t 

e Sl,0, g + e ^0,1, gj J ' 

Re[^] + ilm[^]l pf J, (D7c) 
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and it is possible to construct a reduced linear QTE for Using Eq. (|7.1ip and normalizing, the above QTE gives 
the qubit in the dispersive frame Eq. (|7.9| with measurement record given by Eq. (|7.10p . 
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